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Abstract 

Phase masks coronagraphs can be seen as linear systems 
that spatially redistribute, in the pupil plane, the energy 
collected by the telescope. Most of the on-axis light must 
ideally be rejected outside the aperture so as to be blocked 
with a Lyot stop, while almost all of the off-axis light must 
go through it. The unobstructed circular apertures of off- 
axis telescopes make this possible but all of the major tele- 
scopes are however on-axis and the performance of these 
coronagraphs is dramatically reduced by their central ob- 
structions. 

Their performance can be restored by using an additional 
optimally designed apodizer that changes the amplitude 
in the first pupil plane so that the on-axis light is rejected 
outside the obstructed aperture of the on-axis telescope. 
An apodizer is assumed to be located in a pupil plane, 
a phase mask in a subsequent image plane, and a Lyot 
stop in a reimaged pupil plane. The numerical optimiza- 
tion model is built by maximizing the apodizer's trans- 
mission while setting constraints on the extremum values 
of the electric field that the Lyot stop does not block. The 
coronagraphic image is compared to what a non-apodized 
phase mask coronagraph provides and an analysis is made 
of the trade-offs that exist between the transmission of the 
apodizer and the properties of the Lyot stop. 
The existence of a solution and the transmission of the 
mask depend on the geometries of the aperture and of the 
Lyot stop, and on the constraints that are set on the on- 
axis attenuation. The system's throughput is a concave 
function of the Lyot stop transmission. In the case of 
a VLT-like aperture, optimal apodizers with a transmis- 
sion of 16% to 92% associated with a four-quadrant phase 



mask provide contrast as low as a few 1CT 10 at 1 A/D from 
the star. The system's maximum throughput is about 64%, 
for an apodizer with an 88% transmission and a Lyot stop 
with a 69% transmission. It is showed that optimizing 
apodizers for a vortex phase mask requires computation 
times much longer than in the previous case, and no result 
is presented for this mask. 

It is demonstrated that apodizers can be successfully opti- 
mized to allow phase mask coronagraphs to be used with 
the full aperture of on-axis telescopes while delivering 
contrast as low as or even lower than what they could pro- 
vide by themselves with off-axis telescopes. 

1 Introduction 

Over the past 20 years, indirect detection meth- 
ods have provided us with a now exhaustive cat- 
alogue of confirmed companions (see for example 
http://exoplanet . eu/catalog/) and an even greater 
list of potential companions (for the Kepler can- 
didates, see [http : //archive . stsci . edu/kepl er/| 
planet_candidates.html). One of the most striking 
conclusion so far is that the fraction of rocky planets (in- 
cluding larger and/or heavier Earths) is larger than previ- 
ously anticipated! Bo rucki et al.|20Tl ; Batalha et al.|2012| l. 
It is however only valid for planets with short periods, 
thus at a small distance from their star (about 0.25 AU in 
the case of a solar-type star). 

Several indirect detection methods exist: radial ve- 
locimetry, transits, astrometry, micro-lensing, pulsar tim- 
ing. Transits have been used to probe exoplanets atmo- 
spheres (B urrows et a l. 2006), and the presence of basic 
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molecules such as water, methane, or carbon dioxide has 
been detected in the atmosphere of 5 planets. 

While direct imaging methods cannot currently be used 
to characterize the planets detected by transits, they can 
be used to look for larger planets, further away from their 
star. Some instruments are already being used (Hinkley 
|etal.|201ip . 

The spatial and spectral resolutions of near-future di- 
rect imaging near-infrared instruments will make it pos- 
sible to study planets around nearby stars in much finer 
details ( |Macintosh et al.|2008HBeuzit et al.|2010[|Beich- 
man et al.||2010||Lagage & European MIRI Team||2010[ 
Sivaramakrishna n et al.|2010||Enya & S. W. Group 201 l| l. 

Many instruments have been proposed to detect and an- 
alyze planet light by direct imaging. In this paper, phase 
mask coronagraphs and apodized coronagraphs are pri- 
marily considered. The stellar Lyot coronagraph (Lyot 



1932; Kenkni ghtp977]l was first m odified by Roddier & 
Roddier ( [Roddier & Roddier|1997| ). The focal amplitude 
mask was replaced by a focal phase mask, a disk centered 
on the star that induces a n phase shift. Destructive in- 
terferences in the following pupil plane lead to a partial 
cancellation of the star light. It was shown later that this 
null can be made total if a spheroidal prolate apodizer is 
used in addition ( Aime et al. 2 002| ). 

Rouan proposed another phase mask, the four quadrant 
phase mask (4QPM) to cancel the on-axis starlight en- 



tirely, without using any apodization (Rouan et al. 2000). 
This mask induces a n phase shift of the electric field in 
two opposite quadrants of the image plane, attenuating 
slightly the light coming from the planets if it happens to 
fall on the edges of the quadrants. Eight octants phase 
mask (80PM) have also been studied. While their geom- 
etry reduces the discovery area in the image plane, they 
are also less sensitive to low-order aberrations (Murakami] 
|et al.|2008[|Carlotti et al.|2Q09] > than the 4QPM. 

These phase mask coronagraphs may suffer from chro- 
maticity: n phase shifts are usually induced by laying 
on the surface of a plane-parallel window a thin trans- 
parent film that is then selectively etched ( |Riaud et al.| 
2003). The phase shift that is induced depends on the 



sist in using semi-achromatic polarizing elements (half- 
wave plates) ( |Mawet et al.|2003| |Boccaletti et al.|2008| ), 
or sub-wavelength gratings (Mawet et al. 2005 ). A Mach- 
Zehnder interferometer can also be used to induce an 
achromatic n phase shift ( |Aime et al. 2007 1, although 
many more optical elements are required than in the other 
designs. 

Finally, |Mawet et al.| ( |2005[ ) presented the vector vor- 
tex coronagraph (VVC). The vortex phase mask (VPM) 
induces an azimuthal phase shift, from to 216, where 
is the azimuth, and / the topological charge, an even in- 
teger strictly greater than 0. The VVC also cancels the 
on-axis star light entirely, but contrary to the 4QPM, the 
geometry of the VPM does not reduce the discovery area. 
Several manufacturing techniques are considered for the 
VVC (liquid crystal polymers, sub-wavelength gratings, 
and photonics crystals), and a partial achromatization is 
possible (for more details see Mawet et al. ( 201 1[ |). 

The Roddier & Roddier mask can be modified to im- 
prove its overall performance and its achromaticity by 
adding a phase-shifting annulus to the central mask to 
form a dual-zone mask (Soummer et al. 2003). Like 
the Roddier & Roddier mask, the dual-zone mask can be 
combined with an apodization, which can induce a per- 
fect null in a broadband as long as its transmission varies 
chromatically in a specific manner (N'diaye et al. 2012). 

All phase masks have been extensively studied through 
numerical simulations, and many of them through lab- 



thickness of the film, its refraction index and on the wave- 
length. A first solution to this problem is to use N masks 
designed as many different wavelengths in as many suc- 
cessive image planes (Galicher et al. 201 1) . Other pos- 



sibilities to induce a partial achromatic phase shift con- 



oratory experiments ( |Mawet et al. 2009 Baudoz et al. 
2010). On-sky observations have been conducted with 
the 4QPM coronagraph (Gratadour et al. 2005| ) and the 
VVC ( |Mawet et al.|2010| ), however all three masks have 
been designed for circular clear apertures, not centrally 
obstructed apertures. As a consequence their performance 
has been greatly limited. The VVC has for example been 
tested on a 1 .5m subaperture of the 5m Palomar telescope, 
reducing its effective resolution to a third of its value, and 
its transmission to a ninth. It was nevertheless able to ob- 
serve successfully 3 of the 4 planets of the HR8799 plan- 
etary system ( Serabyn et al.|2010] ) (more information on 
the discovery of these planets can be found in Marois et al. 
( |2008 2010| )). In the specific case of this observation, us- 
ing a subaperture also had the advantage of mitigating the 
impact of wave front errors. 

Phase masks coronagraphs will be used in three major 
high-contrast instruments. On the ground, the spectro- 
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polarimetric high-contrast exoplanet research (SPHERE) 
instrument will use three 4QPM designed for the J, H and 
K bands (Beuzi fet al.|2010"] l. In space, the mid-infrared 
instrument of the James Webb space telescope (JWST- 
MIRI) will also use a series of 4QPM ( |Cavarroc et al.| 



2008 1. SPHERE should see its first light in 2013, while 



the JWST should be launched in 2018. A partially trans- 
missive Roddier & Roddier mask is currently studied for 
the SCExAO instrument at the Subaru telescope ( |Marti 



|nache et al.||201l| >, where it would be used in the focal 
plane of a phase-induced amplitude apodization (PIAA) 
coronagraph ( Guyon|2003[|Guyon et al. 2010|l. Still at the 
Subaru telescope, 4QPM, 8QPM ( [Murakami et al.|20T0] > 
and VPM are also likely to be tested. 

A possible solution to the problem of on-axis observa- 
tions, apart from the use of subapertures, involves two- 



stages optical layouts. For instance, two VPM (Serabyn 



[eTaLl[20TTT l or two 4QPM ( |Galicher et aLl|20TT 



can be 

set in series (in two successive image planes). This re- 
duces the impact of the central obstruction. It also in- 
volves twice as many optics, each of them bringing addi- 
tional aberrations to the wavefront, and making the align- 
ment more difficult than is a single stage layout. 

This paper presents a different solution to on-axis ob- 
servations with phase mask coronagraphs. Apodizers can 
be optimally designed in two dimensions to cancel the on- 
axis light of a star when used with a 4QPM or VPM. The 
apodizer is set in a pupil plane that comes before the im- 
age plane in which the phase mask is located. The opti- 
mization of the apodizer is in many ways similar to the 
optimization of 2 dimensional binary apodizers, i.e. 2D 
shaped pupils. 

While binary pupil masks have long been used by as- 



tronomers (Schultz & Frazier 19831, their optimization 
for high-contrast has however first been proposed about 
ten years ago (Spergel & Kasdin 2001). Several differ- 



ent types of masks have been designed since ( Vanderbei 



et al.|20 03 ; Kas din et al.|2007 1, and they share properties 



with external occulters. While initially optimized in one 
dimension, it was shown more recently that they can be 
optimized in two dimensions (Carl otti et al.||20lT l. This 
makes it possible to tailor the pupil plane transmission 
pixel by pixel. As a consequence apodizers can be com- 
puted for any aperture, including centrally obstructed and 
segmented apertures. 

The optimization problem consists in maximizing the 



transmission of the apodizer under constraints set on the 
extremum values of the amplitude of the electric field in 
chosen regions of the image plane. Because the problem 
is convex and the electric field is computed through a lin- 
ear operation (a Fourier transform), a unique solution ex- 
ists and can be found efficiently. To limit the computation 
time required by these numerical optimizations, the num- 
ber of resolution elements along both axes of the pupil 
plane is usually kept under a thousand, and under a few 
tens in the image plane. 

Computing a shaped pupil that restores the perfor- 
mance of a phase mask coronagraph is also possible, 
though at first glance the complexity of the computations 
appears to be much greater since one or two more Fourier 
transforms must be computed. Constraints are set on the 
extremum values of the amplitude of the electric field, ei- 
ther in the Lyot plane (within the limits of a Lyot stop) or 
in the final image plane. 

The problem can however be simplified by noting that 
the operator that transforms the electric field in the first 
pupil plane into the electric field in the second pupil plane 
is linear. It is thus not mandatory to explicitly compute the 
two successive Fourier transforms. This can drastically 
reduce the complexity of the computation. More impor- 
tantly, it prevents computation errors to spread from one 
plane to another and eliminates the constraint on the sam- 
pling of the electric field in the first image plane. 

Even if the electric field in the intermediate image plane 
is not explicitly computed, the finite angular extension of 
the image plane must be specified. Phase mask corona- 
graphs require the phase shift to be applied on the largest 
possible area of the image plane, usually within a cir- 
cle with a minimum radius of 32 or 64A/D ( |Riaud et al.| 
|200TT >. 

Sec. |2]derives the equations that are later used to com- 
pute the electric field in the Lyot plane as a function of the 
apodizing pattern used in the entrance pupil plane, and 
of the phase-shift induced by the mask in the first pupil 
plane. The case of the 4QPM and of the VPM are con- 
sidered separately. Section [3] gives the details of the dis- 
cretization of the calculations and the construction of the 
optimization problem. Sec. [4] analyzes the impacts of the 
central obstruction and of the finite size of the focal mask 
on the high-contrast performance of a 4QPM and a VPM. 
The results of the optimizations are presented in the same 
section, before being discussed in Sec. [5] 
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Figure 1: Schematic drawing of the optical layout. The 
successive pupil and image planes A, B, C and D are the 
locations of the apodizer, the phase mask, the Lyot stop 
and the camera. The plano-convex lenses were chosen 
to denote the use of a focusing element. They could for 
example be replaced with achromatic doublets, spherical 
or off-axis parabolic mirrors. 

2 Fourier optics formalism 

The optical layout is represented in Fig. [T] Four planes 
are indicated, noted with letters A,B,C and D. They cor- 
respond to the successive pupil and image planes, where 
the apodizer, the phase mask, the Lyot stop and the cam- 
era are respectively located. The telescope's pupil is as- 
sumed to have a diameter D. The observation wavelength 
is noted A, and the letter F denotes the focal length of the 
lenses or mirrors used to reimage the various pupil and 
image planes. 

The apodizer and the telescope aperture are assumed to 
be symmetric with respect to the x and y axes, and since 
their transmission is real, the Fourier transform E(u, v) of 
the apodizer A(x,y) is a real, even function (with the ex- 
ception of a scalar uniform phase factor). Its expression 
is given as follows, in angular units: 



P(x, y) = -OF e 2 ' mF,A f f E(u, v) x M FQ (u, v) e muSt+v ®du dv 
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2.1 Case of the four- quadrant phase mask 

The 4QPM can be described through the odd function 
MfQ(u, v) = sign{u) x sign{v), and the electric field in 
the reimaged pupil plane is the inverse Fourier transform 
of E(u, v) x M FQ (u, v): 



(2) 



For the moment, the dimensions of the focal phase 
mask are assumed to be infinite. After a few algebraic 
transformations, this expression can be reduced to: 



P(x, y) = 4UF e 2inF,x jj E(u, v) sin (2mtx) sin (2nvy) du dv 

(3) 

Plugging in the expression for E(u,v) as a function of 
A(x,y) leads to: 



P(x, y) = 16 e 4i * m J J J J A(x, y) cos (2mix) cos (2ttvv) 

x sin {2nux) sin {2nvy) dx dy du dv 

(4) 

The integrations over u and v can be considered sepa- 
rately: 



P(x,y) = \6e AinFIA I A(;t,;y)[ cos (2nux) sin (2nux) du 
dy. J Jo Jo 

x I cos (2nvy) sin (2nvy) dv dx dy 
Jo 

(5) 

or equivalently: 

rrD/2 

P(x, y) = 16 e 4mF/A I A(x, y) C^x, x) CUy, 50 dx dy 

(6) 
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where the function Cco(a,/3) that has been introduced is 
defined as: 



be represented by the product of a function of u and a 
function of v: 



Coo(a,/3) — I cos (2mia) sin (2mi0) du (7) 
Jo 

While it is not possible to find an analytical expres- 
sion for Cm if the integration domain is semi-infinite, 
it becomes possible for a finite domain. In fact phase 
masks can only cover a finite area of the image plane, 
and this restriction - instead of being an artificial assump- 
tion - is mandatory. It also means that any light outside 
this finite domain is blocked. This loss of throughput is 
however negligible (for example, more than 99% of the 
Airy pattern's energy is concentrated inside a circle of ra- 
dius 20A/D, and more than 90% inside a circle of radius 
2A/D). The expressions of Cl (where L is the angular 
half-width covered by the mask) is written, in units of 
A/D: 



M v ( ff) = cos(0O + i sin(6»0 



= cos'(0)g/*Qtan*(0) 



(9) 



Because the phase induced by the mask depends on 6, the 
integration in plane B is done with respect to the polar 
coordinates p and 9 instead of u and v: 



P(x,y) — iAFe 
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(10) 



Plugging in the expression for the electric field E(u, v) as 
a function of A(x,y) gives: 



Ci(a,y3) = f cos (2nua) sin (2mi0) du 
Jo 

sin 2 (7rL(a + /?)) sin 2 (7rL(a - /?)) 



(8) 
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This concludes our analysis of the mathematical for- 
malism used for the four-quadrant phase mask. Indeed, 
to compute the electric field P{x, y) one only has to plug 
in the expressions of Ci found in Eq(8] Note that in this 
case the mask cover a square region of the image plane 
centered on its origin. 

A closed form for the Fourier transform of P(x, y), and 
therefore the attenuated image of the star, could be found 
as long as no Lyot stop is introduced in that plane. As 
using a Lyot stop appears however mandatory, this addi- 
tional step is not explored further. 

2.2 Case of the vortex phase mask 

The VPM can be represented mathematically through the 
complex function M v (6) = e' eI , where /, the topologi- 
cal charge, is an even integer strictly greater than 0, and 
tan(fl) = v/u. The real part of this function is even (with 
respect to u and v), while its imaginary part is odd. Con- 
trary to the 4QPM, the phase induced by the VPM cannot 
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Where F(x,x,y,y) denotes the double integral with re- 



spect to p and 0. One can notice that Eq. 1 1 can also be 
written as the convolution of A(x,y) and a kernel K(x,y) 
which is the Fourier transform of the mask: 



K(x,y) 



o Jo 



p e im e -2mp(,x cos(e )+ y since)) dpd0 (1 2) 



In any case, the integral with respect to p has a closed 
form if the domain of integration is finite and not semi- 
infinite: 

Jo ^ a 

A closed form can then be sought for the integral with 
respect to 6. We only consider here the case of the topo- 
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logical charge 1 = 2: 



back at Eq|6] this leads to writing that: 
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The electric field P(x,y) can be computed using EqfTT] 
the expression of F(x, x,y,y) being given in Eq{l4] Note 
that contrary to the case of the four-quadrant phase mask, 
and contrary to the case of a classical Fourier transform, 
the integrals with respect to x and y cannot be separated. 
This is a major obstacle to a fast computation of the 
electric field in plane C. The symmetries of the apodiza- 
tion pattern A(x,y) and of the real and imaginary parts 
of F(x, x,y,y) can however be used to integrate over a 
smaller domain of plane A: 



(16) 



Assuming a uniform distribution of points along both the 

p 

2N, ■ 



x and y axes, Ax and Ay are replaced with -Mr . Computing 



the value of the electric field for the points in plane C 
is done with a total of N*N% operations. This number 
can be reduced by proceeding to this integration in two 
successive steps: 
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(17) 



(15) 



While the integration is then done over a fourth of the 
pupil, the function F(x, x,y,y) has to be evaluated four 
times if no simplification can be found and the total num- 
ber of operation remains the same as in the previous case. 
The number of constraints though is reduced by four. 



3 Discretization & Optimization 

To be suitable for computations and numerical optimiza- 
tions, the equations derived previously must be converted 
to a discrete form. 



3.1 Discretization 

The axes in plane A are sampled with Ny points, and with 
N2 points in plane C. In the case of the 4QPM, and looking 



And the complexity becomes N2 + N^N\ . 

In the case of the VPM, using Eq[TT] the integration is 
slightly different: 



D 2 Ni m 

P(x k ,yi) = — I e 4inF/A £ J] Mx t ,yj)F(x i ,x k ,yj,yi) 

i=-N, + l j=-Ni + l 

(18) 



Using Eq[l3] instead of Eq[TT] does not change the com- 
plexity, which equals 2 x 16 x N^Nf in both cases. In- 
deed, F(x, x, y, y) is a complex function, and the real and 
imaginary parts must be considered separately in the op- 
timization, which increases the number of computations 
by a factor two. Furthermore the entire pupil plane must 
be considered here, and four times as many points are re- 
quired. It was previously showed that the parity proper- 
ties of F can be used to compute the field using only one 
quadrant of the pupil, however the actual number of com- 
putations remains the same. Finally, separating the double 
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integral in two successive integrals actually increases the 
complexity to 2 x (16 x N 2 N\ + 8 x N\N\) (in this case 
Ptemp would be a function of 3 variables, and not two). 

One may wonder if it would not be more efficient to 
compute the electric field in plane C through two succes- 
sive Fourier transforms. The total complexity would go 
as low as N\M + N X M 2 + 2 x (M 2 N 2 + MN%), where M 
would be the number of points used to sample the electric 
field along both axes of plane B. This complexity assumes 
a 2-axes symmetry in the entrance pupil plane, and relies 
on the physics of the VPM, so that the electric field can be 
computed in only one quadrant of planes A, B and C. The 
2-axes symmetry is a property shared by many telescopes 
such as the VLT, the Subaru telescope and the Gemini 
telescope. Thanks to the properties of the Fourier trans- 
form, the same 2-axes symmetry is found in the electric 
field in plane B, and the VPM then creates in plane C 
an electric field whose real part presents the same sym- 
metry, and an imaginary part that presents a 2-axes anti- 
symmetry. Hence, it is possible to use these properties 
to compute the electric field in only one quadrant of each 
planes, thus reducing the number of computations by a 
factor 4 in each transform. 

Realistic values for N\ , N2 and M must be chosen to as- 
sess the interest of this method over the previous one. If 
the electric field is not computed explicitly, only the pupil 
planes must be correctly sampled. Both the 4QPM and 
the VPM, used with a circular clear aperture, produce rip- 
ples in the reimaged pupil plane with a periodicity that 
is the invert of the angular extension of the focal phase 
mask. Hence, if the mask extends up to L A/D, 2 L rip- 
ples appear in plane C. Assuming that the patterns that 
the apodizers produce present the same periodicity, the 
value of N\ and N2 must be chosen accordingly, i.e. they 
must be large enough to correctly sample this signal. To 
obey the Nyquist-Shannon sampling theorem, each period 
should be sampled with 2 points, which gives ^=2x1 
(only one quadrant of the pupil plane is actually sampled). 
Thus, observing up to 32 A/D implies that N2 must at least 
equals 64. It might be preferable to consider a higher 
number Ni of points in plane A since, when optimizing 
shaped pupils in two dimensions, the quasi-binary trans- 
missions that are found can be rounded without changing 
the PSF properties as long as N\ is large enough (a few 
hundred points for contrast as low as 10 and a thou- 
sand points for contrast values as low as 1CT 9 ). While 



algorithms exist to turn a gray apodization pattern into a 
binary one - as it is for example done when computing 
microdots apodizers (Martinez et al. 2009]) - a native bi- 
nary solution would have the advantage of preventing any 
potential impact on the performance of the mask that the 
non-binary to binary conversion process would have. 

In the case where N\ — 128 and N2 = 64, the complex- 
ity is 1.6 x 10 6 for the 4QPM and 1.3 x 10 8 for the VPM. 
Choosing N\ = 512 and N 2 = 64 makes it 1.9 x 10 7 for 
the 4QPM and 2.1 x 10 9 for the VPM. 

If the electric field in plane B is explicitly sampled, it is 
important to sample each A/D angular distance unit with 
a high enough number of points, as errors due to a poor 
sampling in plane B will propagate to plane C. Although 
the Nyquist-Shannon sampling theorem requires a mini- 
mum of 2 points per A/D, the specificities of the VVC 
demand a higher number of points. It is also likely that 
the higher the topological charge, the higher the number 
of points since the phase variation across the image plane 
becomes steeper. In fact, in the case of a clear aperture, 
a minimum number of 4 points per A/D is necessary to 
obtain contrast values down to 10~ 8 when the charge is 2. 
When the charge is twice as high, 8 points are necessary 
to recover the same contrast. 

If N\ = 128 and N2 = 64, and the same angular exten- 
sion of 32 A/D is chosen in plane B, then the complexity 
of this computation equals 7.3 x 10 6 if 4 points are used to 
sample one A/D, and it becomes 2.3 x 10 7 if 8 points are 
chosen instead. These numbers are respectively 4.6 and 
14.4 as high as what was required with the 4QPM, but 
also about 300 and 100 less than what would be required 
for the VVC if the electric field in plane C was computed 
directly, using EqjT8| 

3.2 Optimization model 

The optimization of the pupil mask is very similar to the 
optimization of 2-dimensions shaped pupils, as described 
in |Carlotti etaL] ( |2"00"9"l > and | Vanderber] ( |20T2| . Though in 
both cases the problems are convex linear, the constraints 
that are specified here are not computed using Fourier 
transforms, but through the transforms described in Eq(6] 
andQI] 

The problem consists in maximizing the total transmis- 
sion of the mask: 

Z*=i ZyLi A(xi,yj)AxAy under the following constraints 
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set on the values of the transmission of the apodizer, and 
on the amplitude of the electric field in the reimaged pupil 
plane, or plane C: 

< A{xi,y<) < 1 , with fey,} e A A 

(19) 

-l(T e/2 < P{x k ,y{) < 10- c/2 , with {x k ,M € A c 

where A A is the region defined by the telescope's aperture, 
Ac is the region defined by the Lyot stop, and c measures 
the attenuation of the intensity in the reimaged pupil plane 
in a logarithmic scale. With a clear aperture, this attenu- 
ation factor is as small as 1CT 5 in the center of the pupil, 
although this number depends on the size of the phase 
mask, and becomes increasingly smaller as the mask be- 
comes larger. Nonetheless, it is a good starting point for 
the optimization problems that must be solved. 

The optimization model is coded using the AMPL lan- 
guage, and it is solved using the LOQO solver developed 
by Robert Vanderbei (see |Vanderbei| ( |1998] l). 

4 Results 

The implementation of the two linear transforms intro- 
duced in Eq|6] and Eq[TT] makes it possible to study the 
effect of the central obstruction on the contrast in the fi- 
nal image plane D, as well the impact of smaller or larger 
mask size L. In the case of the 4QPM it also allows nu- 
merical optimizations to be performed efficiently. 



The computation of the PSFs showed in Fig(3]assumes 
the presence of a Lyot-stop in plane C. While two dif- 
ferent Lyot stops must be used (one for the clear aper- 
ture, and one for the VLT-like aperture), both have a rel- 
ative transmission of about 90% with respect to the main 
aperture. The 4QPM used with the clear aperture cancels 
the starlight and creates a 10~ 9 contrast floor. The VPM 
performs a little better and creates a 10" 10 contrast floor. 
Theoretically, both the 4QPM and the VPM cancels com- 
pletely the on-axis light and this non-infinitely small value 
can be explained by (a) the finite size of the phase mask 
and (b) the necessary discretization of the pupil. 

The performance of both phase masks with the cen- 
trally obstructed aperture of the VLT is significantly 
worse. Although the maximum intensity of the star is re- 
duced by a factor 10 24 ~ 250, the coronagraphic PSF re- 
mains only at about one to two orders of magnitude below 
the telescope's PSF. 

The contrast curves showed in Fig(3]do not take into 
account the partial cancellation of the planet's light when 
it happens to be too close to the star. Following the defi- 
nition proposed in |Guyon~e t al. (2006), the inner working 
angle (IWA) is defined as the distance for which no more 
than 50% of the planet's light is cancelled. In the case 
of the 4QPM, this IWA is theoretically of about IA/D but 
the planet's light is also attenuated if the planet is located 
on one of the phase transition axes of the mask, in which 
case the edge acts as a phase knife coronagraph canceling 
80% of the light. 



4.1 Impact of the central obstruction 

The transforms displayed in Eq|6]and Eq{TT]are here used 
in two basic cases: one is a circular clear aperture corre- 
sponding to an off-axis telescope, and the other is an on- 
axis VLT-like obstructed aperture (the details of the ge- 



ometry were found in Guerri et al. ( 2008 1). The normal- 
ized intensity that can be observed in the entrance pupil 
and the Lyot plane are displayed in Fig. [2] In both cases 
the linear number of points used to discretize the pupil in 
plane A equals 512, and the focal phase masks extend up 
to 128/l/Z). In plane C, whereas all light is rejected out- 
side the clear aperture, some light is diffracted inside the 
centrally obstructed aperture, reducing the attenuation of 
the on-axis light, as it can be seen in images (d) and (f) of 
Fig.0 



4.2 Impact of the finite angular size of the 
mask 

In practice, whether the 4QPM or the VPM is used, the 
phase mask can always only cover a finite area of the im- 
age plane. In the rigorous analytical demonstration of the 
nulling property of the 4QPM, it is however assumed that 
the mask covers the entire image plane, and it is thus in- 
teresting to study the impact a finite angular size of the 
mask has on the PSF in the final image plane. The an- 
gular size L of the mask is a parameter that must be set 
when computing the electric field in the Lyot plane, and 
Eq|6]can be conveniently used to assess the evolution of 
the contrast in plane D as function of L. 

As it is presented in Fig|4j both the mean and the max- 
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Figure 3: Azimuthal averages of the PSF of (a) the 
4QPM coronagraph and (b) the VVC, for the two aper- 
tures shown in Figj2] The blue line corresponds to the 
case of a clear circular aperture, while the red line corre- 
sponds to the case of the VLT aperture. In both cases, a 
1 Lyot stop with a transmission of about 90% was used to 
try and mitigate diffraction effects due to the finite size of 
L the telescope aperture. 
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Figure 2: Fig. (a) & (b): normalized intensity in plane 
A for a clear and and centrally obstructed apertures. Fig. 
(c) & (d): normalized intensity in plane C for the same 
apertures, after a 4QPM is located in plane B, extending 
up to 128/l/D. Fig. (e) & (f): normalized intensity in 
plane C for the same apertures, after a VPM is located in 
plane B (its angular extension is the same). 





Figure 4: Evolution of the mean contrast (solid black line) 
and the maximum contrast (dashed red line) in plane D of 
a 4QPM coronagraph (left) and a VVC (right) used with a 
circular clear aperture, for a mask radius L up to 320 A/ D 
(for the 4QPM) and up to 64A/D (for the VPM). An 80% 
transmission Lyot stop was used in this case, and both 
pupil planes A and C were discretized over 4096 by 4096 
points (for the 4QPM) and over 5 12 by 5 12 points (for the 
VPM). 
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imum contrasts decrease with higher values of L. In the 
case of the 4QPM, the diffraction effects caused by the fi- 
nite size of the mask lead to a 10~ 5 ' 9 maximum contrast 
for L - 16A/D, a 1CT 8 maximum contrast for L = 48A/D, 
and 10~ 10 maximum contrast for L * 140A/D. The case 
of the VPM is rather similar, although the decrease of the 
contrast value is less monotonous. 

These values were computed assuming a clear circular 
aperture without a central obstruction, and without spi- 
ders. In both pupil planes A and C, the electric field 
was computed over 4096 by 4096 points in the case of 
the 4QPM, and over 512 by 512 points in the case of 
the VPM (this smaller number being chosen because of 
the higher complexity of the computation). The mean 
and maximum contrast values were computed from az- 
imuthally averaged PSFs, over an angular range going 
from to 48A/D for L < 64, and from to 0.75L oth- 
erwise (this was done so as to mitigate diffraction effects 
occurring near the outer edge of the mask). The Lyot stop 
was chosen in this case to have an 80% transmission. A 
more transmissive Lyot stop (90% transmission) leads to 
similar but slower decrease, with a 10~ 10 maximum con- 
trast being reached for L =190 A/ D. 

4.3 Results of the optimizations 

A 3.3 GHz processor out of a 24 cores machine, with a 
total RAM of 192 GB was used to solve the optimiza- 
tion problems. Since in practice this machine is shared 
with many users, only up to 8GB are available for one 
optimization process. Because of its lower complexity, 
the case of the 4QPM was considered first. Five differ- 
ent apodizers have been computed for this phase mask 
(see Fig|5]) and for a VLT-like aperture. In the case of 
the VPM, the complexity of the optimization process is 
not much higher, but too much RAM was required and no 
results have yet been obtained for this mask. 

Every apodizer has been designed to attenuate the 
diffracted light inside the area defined by a specific Lyot 
stop (see Fig|6]). The five Lyot stops were computed by 
oversizing the central obstruction and the spiders of the 
telescope, while undersizing its outer diameter. To do 
that, an image filtering process was used: an array de- 
scribing the telescope aperture is convolved with unit ma- 
trix. The result is a blurred image of the aperture. Round- 
ing it creates a Lyot stop. The transmission and resolution 



losses are adjusted by the relative size of the unit matrix 
and the aperture. An array of 1024 by 1024 points was 
used for the aperture to accurately change the transmis- 
sion in the pupil plane. The Lyot stops were then interpo- 
lated in two dimensions to reduce their size. The transmis- 
sions Tis of the five Lyot stops range from 59% to 83%. 
These transmissions are defined as the ratio of the total 
intensity going through the masks with the total intensity 
going through the initial aperture of the telescope. 

Indeed, each apodizer is optimized over an array of 128 
by 128 pixels covering one quadrant of plane A. In ev- 
ery cases, the angular extension L of the phase mask is 
32 A/D. The electric field in plane C is also computed 
over one quadrant, and with the same number of points 
as in plane A. An average time of a day and a half were 
necessary to compute each apodizer. 

The effect of the apodized 4QPM on the spatial distri- 
bution of the energy density in the Lyot plane is showed in 
Fig|7]for mask #4. On-axis light is concentrated outside 
the transmissive area of the Lyot stop and is thus almost 
entirely blocked by it. On the contrary, in the unapodized 
case, light is spread in these regions and the extinction of 
the star is rather poor. 

Radial averaged profiles of the PSFs are displayed in 
FigjH] for the five apodizers (two-dimensional contrast 
plots were averaged azimuthally so as to be represented in 
one dimension only). A contrast floor lower than 10~ 8 and 
higher than 10~ 10 is obtained with each mask, the lowest 
ones corresponding to the least transmissive Lyot stops. 
For comparison, the PSF of the unapodized 4QPM is also 
showed. A gain of four to five orders of magnitude is pro- 
vided by using an apodizer. 

A comparison of these results with those displayed in 
Figj4] makes it clear that the apodizer does not only help 
restoring the contrast performance of the 4QPM. While 
a maximum contrast of 10~ 7 should have been obtained 
for mask size L of 32 A/D, maximum contrast as low as 
10 are seen in the PSF showed in Fig §(the maximum 
is looked for over an angular range that extends between 
and 24 A/D, as done it sec |4.2| i. Hence, in addition to 
limiting the impact of the central obstruction, the apodizer 
helps mitigate the diffraction effects due to the finite size 
of the phase mask. 

It is important to notice that these PSFs were computed 
using the 'raw' transmission of the apodizers. While 
the transmissions of the masks tend to be binary, a 5 to 
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7% Rhi„ fraction of the pupil area exhibits a transmission 
greater than 1% and smaller than 99% (see Tab[T]i. FigfTT] 
compares the result of two optimizations that differ by the 
size of the arrays used to discretize the apodizer in plane 
A: the first one has 128 by 128 points while the other has 
256 by 256 points. The apodizers' transmissions are sim- 
ilar, but there are more than 3 times less points with a 
non-binary transmission in the larger array (Rbm * 1.8%). 
This is not unexpected as it has already been noticed in 
2D shaped pupil optimization problems. 

Artificially rounding the transmission of the apodizers 
has a significant impact on the contrast. In the case of 
mask #4, as displayed in Fig[TT] the contrast floor is be- 
tween 10 7 and 10~ 6 . However, the apodizer computed 
for the same Lyot stop transmission, but over four times 
as many points, gives a lower contrast that remains be- 
tween 1CT 8 and 10~ 7 beyond 3A/D, and that does not ex- 
ceed 1CT 6 6 below this angular distance. 

As it is summarized in Tab[T] the respective transmis- 
sions T A of the five apodizers range from 92% to 16%. 
These transmissions are defined in the same way as the 
Lyot stops transmissions were previously defined, as the 
ratio of the total intensity going through the apodizers 
with the total intensity going through the initial aperture 
of the telescope. 

The maximum throughput T max of the system ranges 
from 64% to 5%. T max is not computed as the product of 
Tls an d Ta, as it is the case with most coronagraphs. To 
do so, one must assume that the focal plane mask does not 
have a significant impact on the planet light transmission, 
which fails to be true in the case of the 4QPM. Instead 
T max is measured by injecting an off-axis source in the in- 
strument and measuring its attenuation in the final image 
plane. 

Because of the phase mask, the throughput T(u, v) 
varies across the image plane. Its maximum, T max , 
is found by moving the planet at equal distance from 
the 4QPM phase transitions, hence along a diagonal. 
Fig|9]shows how the normalized throughput T(u, v)/T max 
changes along the diagonals of the masks, and in general 
in the image plane (but only in the specific case of mask 
#4). 



Mask# 


T LS 


T A 


T 

1 max 


IWA 


Rbin 


1 


83% 


16% 


5% 


1.37 A/D 


6% 


2 


79% 


42% 


23% 


1.10 AID 


7% 


3 


74% 


66% 


44% 


1.06 AID 


7% 


4 


69% 


88% 


64% 


1.01 A/D 


6% 


5 


59% 


92% 


56% 


0.99 A/D 


5% 



As it is illustrated in Fig 10 to each Lyot stop corre- 
sponds a different value of T max , which goes through a 
maximum of 64% with mask #4 (Tls = 69%). This max- 
imum value can be roughly estimated by extrapolating 



Table 1 : Main parameters of the five different apodizers 
optimized for a 4QPM coronagraph. Each mask corre- 
sponds to a specific Lyot stop, with a transmission Tls- 
Ta is the transmission of the apodizer. This results in the 
system's total throughput T. The IWA of the coronagraph 
is given in units of A/D. Finally, R/, in is the ratio between 
the non-binary surface of the mask and the total surface 
of mask. The non-binary elements are defined as having 
a transmission greater than 1% and less than 99%. 



from a linear fit based on the evolution of T max with Tls ■ 
The line that can be drawn crosses the line T max = Tls 
for Tls ~ 68%. Since the system's maximum throughput 
T max cannot be larger than the Lyot stop's transmission 
Tls > T = 68% can be thought as an upper-bound estimate 
of the maximum value of T max for this specific telescope 
architecture, and this family of Lyot stops. Still in FigfTU] 
the product of Tls an d Ta is reproduced, as a function of 
T L s ■ ^ can be noticed that this product overestimates the 
value of T max for masks #1, 2 and 3 and underestimates it 
for masks #4 and 5. 

The system's maximum throughput T max is however not 
the only metric one should refer to to evaluate the sys- 
tem's performance. One must also consider the angu- 
lar resolution of the instrument. Using the information 
displayed in Fig. [9] and following [Guyon et al.| ( |2006] l, 
the effective IWAs range between 0.99 and 1.37 A/D (see 
Tab{T]i. The smallest IWA corresponds to the 59% trans- 
mission Lyot stop, and the largest IWA to the 83% trans- 
mission Lyot stop. 

The time required in the case of the larger array equals 
1 18h, which is 3.3 times longer than the time required for 
the smaller array. This is in accordance with our estimate 
of the complexity, as the ratio of the complexities is 3 in 
this case (Ni = 256, N 2 = 128). 
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Figure 5: Transmitted intensities of the five apodizers op- 
timized for the 4QPM at the VLT and a Lyot stop with a 
transmission of (a) 83%, (b) 79%, (c) 74%, (d) 69%, (e) 
59%. White denotes a unit transmission, and black a zero 
transmission. 
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Figure 6: Contours of the VLT aperture (red) and of the 
five Lyot stops (black) with the respective transmissions 
of 83, 79, 74, 69 and 59%. 



5 Discussion & Conclusion 

Both the original four-quadrant phase mask coronagraph 
and the vector vortex coronagraph see their performance 
reduced with centrally obstructed apertures such as the 
VLT's, and the contrast they provide is only 1 to 2 orders 
of magnitude below the PSF of the telescope. On the con- 
trary, when used with clear circular apertures, these coro- 
nagraphs could provide contrast as low as 10" 10 at only 
IA/D from the star. 

Solutions have been proposed in the past that consist 
in either using a circular subaperture, or using a two- 
stages optical layout. The first solution leads to an im- 
portant resolution loss and an even higher transmission 
loss. The second solution offers a more satisfying result, 
but requires twice as many optical components, which can 
limit the system's performance because of alignment is- 
sues, additional wavefront aberrations and chromaticity 
effects. 

It was shown that optimal apodizers can help restore 
the performance of these coronagraphs when used with an 
on-axis telescope, and that the diffraction effects caused 
by the finite size of the focal phase mask are also miti- 
gated thanks to the apodizer. In addition to the already 
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Figure 7: Log scale intensities in the Lyot plane (a) with- 
out the Lyot stop and without apodization, (b) with the 
Lyot stop and without apodization, (c) without the Lyot 
stop and with apodization, and (d) with the Lyot stop and 
apodization. Mask #4 is used here. It is designed for an 
69% transmission Lyot stop. 



Figure 8: Azimuthally averaged contrast obtained with 
the apodizers designed a Lyot stop with a transmission of 
(a) 83%, (b) 79%, (c) 74%, (d) 69%. In each figure the 
blue line denotes the PSF of the telescope (including the 
presence of the Lyot stop), the black line corresponds to 
the PSF of the 4QPM without apodization, and the red 
line corresponds to the PSF of the apodized 4QPM. 
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Figure 9: (a): Normalized transmission of the apodized 
4QPM for the five different Lyot stops, and for a planet 
located, along the diagonal, between and 8 A/D. Sorted 
by increasing transmission, the curves correspond to the 
83%, 79%, 74%, 59% and 69% transmission Lyot stops, 
(b): two-dimensional map of the normalized transmission 
in the case of mask #4 (Lyot stop transmission: 69%). 



Apodizalion - LS: 69% - N=128 Apodizalion - LS: 69% - N=256 




(a) (b) 



Figure 1 1 : Transmitted intensities of two apodizers opti- 
mized over arrays of different size, for the 4QPM at the 
VLT and an 69% transmission Lyot stop. Only one quad- 
rant of the pupil plane is showed. The number of points 
in these arrays is (a) 128 by 128 and (b) 256 by 256. 
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Figure 10: Evolution of the system's throughput T max as 
a function of the Lyot stop transmission Tls ■ The vertical 
and horizontal dashed lines indicate the coordinates of the 
intersection point between the L = Lls line and the linear 
trend of the four points of the solid curve plotted at the 
right side of the figure. The dashed curve represents the 
evolution of the product Tls x Ta as function of Tls ■ 



Figure 12: Azimuthally averaged contrast plots obtained 
with the apodizer designed for a 69% transmission Lyot 
stop. As in Fig j8] the blue line denotes the PSF of the 
telescope (including the presence of the Lyot stop) while 
the black line corresponds to the PSF of the 4QPM with- 
out apodization. The three red lines correspond to the PSF 
of the apodized 4QPM. The solid line is used for the PSF 
of the 'raw' mask (i.e. non-binary), and the dashed and 
dotted lines are used respectively for the PSFs of the bi- 
nary masks computed over 256 by 256 points and 128 by 
128 points. 
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existing phase mask and Lyot stop, an apodizer is located 
in an upstream pupil plane. No other additional optical 
components are required. Any aperture can be used as a 
basis for an optimization, although symmetric apertures 
have the advantage of reducing the number of points and 
thus the complexity of the computations. 

The numerical optimization of the apodizer's transmis- 
sion can be done in two dimensions as it is already done 
for shaped pupils. The complexity of the problem is how- 
ever different whether a 4QPM or a VPM is used. Dif- 
ferent cases require different computations methods. In 
the case of the 4QPM, the direct computation of the elec- 
tric field in the Lyot plane appears to be a better choice. 
Indeed, sampling errors in plane B are avoided since the 
electric field need not be computed explicitly in this plane. 
In the case of the VPM however, the computation of the 
same electric field should better be done through two suc- 
cessive Fourier transforms (from pupil plane A to im- 
age plane B, and then to pupil plane C). While 1 .5 days 
are necessary to compute an apodizer over a 128 by 128 
matrix for the 4QPM, the same optimization would take 
about 3 to 9 times longer for the VPM, depending on the 
sampling resolution of plane B. 

The apodizers that are presented here are designed only 
for the 4QPM. A VLT-like aperture was chosen and dif- 
ferent Lyot stop sizes were considered. For each Lyot- 
stop, an apodizer is found by solving a linear convex opti- 
mization problem. The total throughput of the system is a 
function of the Lyot stop's transmission. It goes through 
a maximum for which an upper-bound of about 68% is 
found for a 68% transmission Lyot stop (and in the partic- 
ular case of the VLT). In practice, i.e. when optimizing an 
apodizer for this specific Lyot stop, the system's through- 
put that is found is 64% and is given by an apodizer with 
an 89% transmission, associated to a 69% transmission 
Lyot stop. 

As observed previously with fully optimized shaped 
pupils, the masks that have been computed tend to have 
a binary transmission: only 5 to 7% of the pupil sur- 
face has a transmission higher than 1% and smaller than 
99%. This number goes down to 1.7% when the size of 
the array is doubled, from 128 by 128 to 256 by 256. 
It is very likely that masks optimized over even larger 
arrays would have an even smaller non-binary transmis- 
sion ratio. Computing a mask over a 512 by 512 array 
would take an estimated time of 15 days. In the case of 



full two-dimensional shaped pupil optimization, array's 
length of half a thousand points have proven to be large 
enough to obtain quasi-binary apodizers whose transmis- 
sion, w hen artificially roun ded, provides contrast down 
to 10~ 9 ( |Carlotti et al.|201 1 1. Contrast values are signifi- 
cantly increased when artificially rounding the transmis- 
sion of the apodizers: up to about 10~ 6 for a 128 by 128 
mask, and up to about 10~ 7 for a 256 by 256 mask instead 
of a few 10 for mask #4. Optimizing masks over larger 
arrays is mandatory, and this point will be addressed in 
the near future, thanks to a new computer, dedicated to 
weeks-long computations. 

Binary masks can efficiently be manufactured using 
photolithography: a thin metallic coating can be deposited 
on a glass substrate and then selectively etched. The sub- 
strate adds chromaticity and wave front aberrations. If the 
mask is structurally connected, it can be etched out of 
a thin silicon wafer coated with metal, eliminating these 
disadvantages. 

Gray apodizers have also been produced, either by de- 
positing on a substrate a metallic layer with a spatially 
varying thickness, or by using high energy beam sensitive 
(HEBS) glass. These apodizers are however considered 
too chromatic (Sivaramakrish nan et al.|2 009). 

Apart from mask #1, every mask computed for the 
VLT-like aperture would require the use of a substrate. 
While a free-standing apodizer would certainly have its 
advantages, the apodized coronagraphs of both SPHERE 
and GPI use substrate-based apodizers. As presented in 
Eny a et aL] (|2008 ), substrate-based shaped-pupils with 
anti-reflection coatings on both sides have achieved 8 x 
10" 8 contrast in the visible. SPHERE and GPI both work 
in the near-infrared where, given the same surface aberra- 
tions, a substrate would create weaker wavefront aberra- 
tions than in the visible, which could potentially lead to 
even lower contrast values. 

The possibility also exists that MEMS or MOEMS de- 
vices could be used as adaptive binary apodizers, either 
in transmission (using a technology close to the one used 
in JWST's micro-shutter arrays) or in reflection (such as 
the micro-mirror arrays used in many video-projectors). 
Apart from offering substrate-free apodizers, it would also 
allow a direct amplitude control of the wavefront. 

Apodizers optimized for a 4QPM could potentially 
complement SPHERE'S 4QPM coronagraphs, increasing 
the achievable contrast while offering the same very small 



15 



resolution. 

No apodizer could yet be optimized for the vector vor- 
tex coronagraph. The VVC has however a substantial ad- 
vantage over the 4QPM since it does not reduce the dis- 
covery area. Both the charge 2 VVC and the 4QPM are 
very sensitive to finite stellar size and jitter. This sensi- 
tivity can be reduced if a charge 4 VVC or an 8 octants 
mask is used instead, although the 8 octants reduces the 
discovery area even more than the 4QPM does. 

There is no inherent reason that makes the optimiza- 
tion of apodizers for the VVC impossible. Shortening the 
computation time would however greatly help this pro- 
cess. To that end it will be necessary to either simplify the 
problem further, possibly by assuming circular-symmetric 
apertures and apodizers, or to solve the optimization prob- 
lem in a more efficient manner, for example by taking ad- 
vantage of the multiple cores of a computer. 

Finally, the joint optimization of the apodizer and the 
Lyot stop is also a possibility that will be explored. 
Though more computationally intensive, it would de facto 
maximize the transmission of the system. 
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